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Abstract 

We review recent work [U El [3] on the Hawking radiation of rotating brane black 
holes, as may be produced at the LHC. We outline the methodology for calculating 
the fluxes of particles, energy and angular momentum by spin-0, spin-1/2 and spin- 
1 quantum fields on the brane. We briefly review some of the key features of the 
emission, in particular the changes in the spectra as the number of extra dimensions 
or the angular velocity of the black hole increases. These quantities will be useful 
for accurate simulations of black hole events at the LHC. 



1 Introduction 



The formulation, nearly ten years ago, of theories with large extra dimensions [H [5] has 
led to a highly active programme of research in gravitational physics (see also [6] for some 
early work in this direction). These theories were motivated by the hierarchy problem 
in particle physics. In order to resolve this problem, the fields of the standard model of 
particle physics are constrained to live on a (3 + l)-dimensional brane, which is embedded 
into a higher dimensional bulk (see figure [1]), so that the total number of space-time 
dimensions in the bulk is 4 + n. Only gravitational degrees of freedom are allowed to 
propagate in the bulk. 

*Talk given at the "Workshop on Dynamics and Thermodynamics of Black Holes and Naked Singu- 
larities H", Milan, 10-12 May 2007. 
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Figure 1: Sketch of the brane- world scenario: our universe is a brane, embedded in a 
higher-dimensional bulk. The standard model fields of particle physics are constrained to 
live on the brane; only gravity can propagate in the bulk. 

One consequence of these higher-dimensional theories is that the scale of quantum 
gravity may be very much lower than the Planck scale Mp. For large extra space-like 
dimensions with a compactification radius L, the fundamental scale of quantum gravity in 
the bulk, M^^ is related to the Planck scale (which is now an 'eflFective' quantum gravity 
scale as seen on the brane) as follows [4]: 

Ml = SttM^^+^L^, (1) 

where n is the number of extra dimensions. It is possible, via equation ([T]), to have 
Mp ^ 10^^ GeV and ^ 1 TeV for various combinations of n and L: for example, 
if n = 2 and L ^ 1 mm. This lowering of the bulk scale of quantum gravity by many 
orders of magnitude may be expected to have significant phenomenological implications. 
In particular, particles with energy above M^ will be able to probe quantum gravita- 
tional physics and collisions of such particles (for example, at the LHC), may result in 
strong gravitational phenomena. In this paper we focus on the formation and subsequent 
evaporation of mini black holes in such trans-Planckian collisions [7J, as may occur in 
ground-based colliders [8j or ultra-high energy cosmic ray interactions in the atmosphere 
of the Earth [9\. We do not attempt to survey here the huge literature on these sub- 
jects, but simply refer the reader to the reviews [TOl Ull [12] for detailed discussions and 
comprehensive lists of references. 

The outline of this article is as follows. We begin in section [2] by very briefiy 
reviewing the key processes in the formation and evaporation of brane black holes at 
the LHC. Our focus is one part of the evaporation, the 'spin-down' phase (which will be 
defined in the next section). Our study involves quantum field theory in curved space, 
and the methodology is outlined in section [31 The key features of the results of this 
study, namely the fiuxes of particles, energy and angular momentum for particles with 
spin-0, 1 and 1/2, are discussed in section [H before we present our conclusions in section[5l 
Throughout this paper, the metric has signature (+,—,—,—) and we use units in which 
G = c = h = kB = l' 
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2 Formation and Evaporation of Brane Black Holes 



The formation of a mini-black hole in a trans-Planckian collision at the LHC may be 
understood heuristically as follows (see sketch in figure [2j). Suppose we have two colliding 



Figure 2: Sketch of the formation of a black hole. If two particles collide in such a way 
that their impact factor b is smaller than the Schwarzschild radius corresponding to their 
centre-of-mass energy, then it is expected that a black hole will be formed. 

particles with centre-of-mass energy E. If their impact parameter b is smaller than the 
Schwarzschild radius of a black hole having mass then, classically, one would expect a 
black hole to be formed [I3]. If the bulk scale of quantum gravity is as low as a few TeV, 
then the Schwarzschild radius of a higher-dimensional black hole is of order <C 
which is very much larger than would be anticipated in general relativity. This means 
that it is feasible that colliding particles at the LHC may have impact parameters smaller 
than the corresponding Schwarzschild radius. The black holes formed by this process will 
be higher-dimensional, extending into the bulk. Of course, this simple argument does 
not accurately portray the whole process, which has been extensively studied [71 [HI [TO] . 
Various different calculations have been done of the rate of black hole production at the 
LHC, with controversial results, but including the exciting possibility of many black holes 
per day being formed ^ i8j. 

In this article we consider black holes with masses many times greater than the bulk 
quantum-gravity scale M*. Such black holes can be considered as semi-classical objects, 
and quantum gravity effects can be effectively ignored. The main observable feature of the 
subsequent evolution of the black hole will be its Hawking radiation [H]. This proceeds 
in four stages (the first three of which are sketched in figure [3]). 

1. When the black hole forms, it will be a complicated object, highly asymmetric and 
endowed with non-abelian gauge field hair from the incoming particles. The first 
stage of its evolution is the ^balding ^ phase^ during which the gauge field hair is shed 
and the asymmetries removed by gravitational radiation. It is expected that this 
phase of the evolution will be extremely rapid. 

2. At the end of the 'balding' phase, the black hole has no hair and will be axisym- 
metric, rotating extremely rapidly. During the ^spin-down^ phase^ the black hole 
then evaporates through Hawking [13] and Unruh-Starobinskii [15] radiation, losing 
mostly angular momentum as well as mass. 

3. At the end of the 'spin-down' phase, the black hole has shed all its angular momen- 
tum and is spherically symmetric. Since the black hole is Schwarzschild-like, the 
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Figure 3: Sketch of three stages of the evaporation of the black hole: (a) the 'balding' 
phase, where the black hole loses its asymmetries and any gauge field hair; (b) the 'spin- 
down' phase, where the black hole sheds it angular momentum and some mass; (c) the 
'Schwarzschild' phase, where the black hole is spherically symmetric and loses mass. The 
last unknown stage, the Tlanck' phase, is not sketched. 

next phase is the ^SchwarzschiW phase^ when the black hole continues to evaporate 
via Hawking radiation, but now in a spherical manner. 

4. The black hole can continue to be described semi-classically until its mass is of the 
order of the (higher-dimensional) quantum gravity scale, in which case unknown 
quantum-gravity effects become important. It is thought that this final stage, the 
^ Planck^ phase will be short-lived compared to the 'spin-down' and 'Schwarzschild' 
phases. 

Of these four stages, the 'Schwarzschild' phase is by far the simplest to analyze and 
correspondingly the most work exists in the literature on this part of the evolution [16]. 
Both analytic and numerical studies have been performed, for standard model fields on 
the brane, and scalar and graviton fields both on the brane and in the bulk. References 
[16j also include work on different types of spherically symmetric black hole background 
(including black hole charge, a cosmological constant, and higher-derivative gravity) and 
on the emission of massive as well as massless quantum fields. 

Our focus in the present article is the 'spin-down' phase. After the extensive work 
on the 'Schwarzschild' phase, the 'spin-down' phase has recently been the focus of much 
attention in the literature [H El [3l [17] . In this article we concentrate on the emission of 
quantum fields of spin 0, 1/2 and 1 on the brane, although there is work in literature [IT] 
which includes emission into the bulk. 
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3 Quantum Fields on a Rotating Brane Black Hole 



If we assume that the coUiding particles which form the black hole are constrained to lie 
on an infinitely thin brane, so that they have a non-zero impact parameter only in the 
brane, then the black hole thus produced will have its plane of rotation in the brane. 
During the spin-down phase, the higher-dimensional black hole geometry will therefore 
be of the Myers-Perry form, with a single non-zero angular momentum parameter a: 

-(^ + a^ + ^^^^] sin^ e dip-" - cos^ 9 dn^, (2) 

\^ TiT^ ^ J 

where 

A = r2 2__^ S = r2 + a2cos2^, (3) 

and d^n is the line-element on a unit n-sphere. In this paper we are interested in the 
emission of particles on the brane, and so we first need to determine the eflFective, four- 
dimensional, black hole metric on which the quantum fields propagate. We do this by 
fixing the values of the additional angular co-ordinates used in ([2]) to describe the extra 
dimensions. Effectively we are taking a 'slice' through the higher-dimensional black hole 
geometry, as sketched in figure [H The four-dimensional metric resulting from this process 




Figure 4: Sketch of the geometry of the higher-dimensional black hole. The black hole 
extends into the bulk, and we consider the 'slice' of the black hole which intersects with 
the brane (our universe). 

is: 

ds^ = (i-^)df + ^-^^^dtd^-^dr^-j:de^ 

(2.2. aVsin^^\ . 2/, > 2 /.N 

The location of the black hole event horizon is given by solving the equation A = 0, 
which, for n > 1, gives the unique solution 
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where a* = a/vh. For n = 0, 1, there is a maximum possible value of a which gives a 
real solution to the equation A = 0, but for n > 1 there is no such a priori bound on 
a. However, there is an upper bound on a from the requirement that the black hole be 
formed by the collision of two particles, as shown in figure O namely [l^j 



^(^ + 2). 



(6) 



We wish to study quantum fields of spin 0, 1/2 and 1 on the background (|4j). To do 
this, we employ the Newman- Penrose formalism [19], and follow the approach of Teukolsky 
[20]. As in Teukolsky's analysis of the perturbations of Kerr black holes, here we find that 
the wave equations can be written as one single 'master' equation, with the 'helicity' 
h = (+5, —s) of the field as a parameter in the equation [TTl [20]: 
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(7) 



where = ^hit^T^^d^cj)) represents the spin-field perturbation and the source term. 
The reader is referred to [HI [20] for the details of this derivation and the precise definitions 
of and Th for each spin. The 'master' equation ([7|) is separable for any value of the 
helicity /i, so that its solution can be written as a sum over the Fourier modes: 
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(8) 



where /^aA are the Fourier coefficients and the set of 'quantum' numbers is given by 
A = {£muj}. The radial and angular ODEs resulting from the separation of variables of 
the 'master' equation ([7|) are, respectively: 







and 



d 



dr \ 



( A^+i 



dr 



+ 



d_ 

dx 



- ihKA\r) 



+ Uhur + h{A"{r) - 2)5h\h\ - /^A 



(9) 



(1 - x^)-^ ) + {aujf{x^ - 1) - 2/ia6JX 



1 



x"- 



dx 



+ h^K + "imauu + h 



hSA(x)] 



(10) 



where x = cos 9 and h^A is the constant of separation between the angular and radial 
equations. 

The angular equation ( ITOl ) has solutions which are the spin-weighted spheroidal har- 
monics [21]. Properties of these special functions and numerical methods for finding them 
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can be found in the literature ^ El [3l EI]. The spin- weighted spheroidal harmonics 
are normalized according to 

/ d9 sin 9 hSl{co8 9) = 1. (11) 

The angular equation (fTOl) has to be solved first because it is necessary to determine the 
separation constant h^A which also appears in the radial equation ([9]). Having found the 
separation constant, we then numerically integrate the radial equation ([9]). There are 
some subtleties about how this is done in practice, depending on the dimensionality of 
space-time and the spin of the quantum field. We do not discuss these issues here; further 
details of our numerical methods can be found in [J, El Ej. We solve the radial equation 
([9]) with the boundary conditions that the solutions must be ingoing at the future event 
horizon (this corresponds to the definition of the usual 'in' modes). From our solution 
of the radial equation, a key quantity, which will be essential later in our analysis, is the 
transmission coefficient Ta, which is defined as the ratio of energy fiuxes: 

dE^l'^'^/dt ^ ^ 

with dE^^^^^ /dt and dE^^^^^ /dt denoting, respectively, the incident energy ffux on the black 
hole and the energy ffux transmitted down the event horizon, for a particular ffeld mode 
with quantum numbers A. 

We are interested in the ffuxes of particles, energy and angular momentum from the 
evaporating black hole due to Hawking radiation, the Hawking temperature of the black 
hole ([4j) being 

_ (n + 1) + (n - l)al 

4.(1 + a^)r. • (''^ 

The evaporating black hole is described by the 'past' Unruh vacuum \U~) [22] on the 
metric (|4j). The ffuxes of energy E and angular momentum J can be expressed in terms 
of components of the renormalized stress-energy tensor: 

dE C 

— = / {U- \T'^U-), en sine ded^; 

— = / {U-\T:^\U-),enr'smeded^; (14) 

while the particle ffuxes are given by the expectation values of the appropriate particle 
number operator. Fortunately the stress-energy tensor components in (fT^ do not require 
renormalization (for spin 0, the argument is presented in [22j, and can be generalized to 
higher spins, using the renormalization counterterms presented in [23j). 

In order to construct the state we start by deffning the usual 'in' and 'up' 

modes for the quantum ffeld [22]. We then take a basis of ffeld modes for which the 
'in' modes have positive frequency with respect to co-ordinate time t near inffnity, while 
the 'up' modes have positive frequency with respect to Kruskal time near the past event 
horizon. We do not describe the details here; instead we refer the reader to the relevant 
literature for the cases of spin [H ESI El], spin 1/2 [21 EH] and spin 1 [3, 26j. After 
computing the expectation value of the stress-energy tensor in the state the total 
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fluxes of particles, energy and angular momentum emitted by the black hole per unit time 
and frequency, for each species, are given by [HEllS]: 



cPN 


(1 + Ss,l) 


dtduj 
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(1 + Ss,i) 


dtdcu 


27r 



oo i 

,_E^;?TOT-- (-) 

where 

UJ = uu — mflHi (18) 
and Qh is the angular velocity of the event horizon, given by 

In equations f fTSUTTI ). the positive sign in the denominator is used for fermions, while the 
negative sign gives the correct Planck factor for bosons. Note that there is an additional 
factor of 2 in the gauge boson (spin-1) case, which results from a sum over the parity of 
the field modes. 

The emission rates f fTSUTTI ) describe the different fiuxes emitted by the black hole 
over the whole solid angle Ql = 4:71^ and follow by performing an integration over the 
angular coordinates 9 and (f. If we go one step backwards, we may derive the angular 
distribution of the emitted radiation by displaying the exact dependence of the differential 
rates on the latitudinal angle ^, and write 

d{cose)dtdu; - Att exp (Co/Tn) ± 1 ^ ^-'^^ + ' ^ 

dPE (1 + 5 i) ^ UJ 

d{cose)dtdu = ~^f~ ^ exp(c:;/rH)±l'^^ ^-'^^ + • ^^^^ 

Again, there is an additional factor of 2 in the spin-1 case, arising from a sum over the 
parity of the field modes. Note that in the spin-0 (scalar) case, the two helicities —h 
are equal to zero, so we simply have an additional factor of two. 

We will now outline the key properties of the fiuxes (11511211) for each of the particle 
species. 
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4 Fluxes of Particles, Energy and 
Angular Momentum 



We now briefly outline some of the key features of our results. For further details, and 
a more wide-ranging discussion (including, for example, the total emissions), we refer 
the reader to the literature ^ El [3l |T7]. Detailed studies of the Hawking radiation from 
four-dimensional black holes can be found in [27\ . 



spin spin 1 
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Figure 5: Transmission coeflicient Ta for spin-0 and spin-1 flelds, for the i = 1 = m mode 
with n = 1, for various values of a*. 

The flrst quantity of interest is the transmission coeflicient Ta ( 1121) . This is plotted 
for the mode £ = 1 = m and n = 1 in flgure O for the cases of spin-0 and spin-1 and 
different values of a^. The interesting feature of the plots in flgure [5] is that in both cases 
the transmission coeflicient is negative for small values of cuTh- This eflFect is small for the 
spin-0 mode but signiflcant for the spin-1 mode. This is, of course, the famous super- 
radiance effect: an incoming wave is reflected back to inflnity with amplitude greater than 
it had originally. This effect is observed only for bosons, and does not occur for spin-1/2 
modes. 

We now turn to the observable fluxes of particles, energy and angular momentum. 
Firstly, the particle flux (1T51) is shown in flgures [6] and [3 as a function of cjr/^, for flxed 
n = 1 and varying a* (flgure [6]) and flxed = 1 and varying n (flgure [71). 

There are number of features that can be seen from flgures [6] and [3 

1. As the number of extra dimensions increases, so do the fluxes, mostly because the 
Hawking temperature (1T31) of the black hole increases linearly with n for flxed a*. 

2. For flxed n, as a* increases, the spectra become broader and flatter. 

3. For smaller values of n, there are oscillations in the particle spectrum, which cor- 
respond to various fleld modes becoming dominant at certain frequencies (a sort of 
resonance effect). 

4. Emission at higher frequencies becomes more important as n or increase. 
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Figure 6: Particle fluxes (fT5l) for spin-O, spin-1/2 and spin-l flelds as a function of uurh^ 
for flxed n = 1 and varying a*. 



spin spin 1/2 




Figure 7: Particle fluxes (fT5l) for spin-O, spin-1/2 and spin-l flelds as a function of ujThj 
for flxed a* = 1 and varying n. 
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5. The peaks for spin-0 and spin-1/2 particles are of the same order of magnitude, but 
the peak in the emission for spin-1 particles is roughly ten times bigger. 

We have also studied how the particle flux depends on the latitudinal angle 9 fl20l). 
and the results can be seen in flgureO In flgure[8l we have cos^ = ±1 along the rotation 





Figure 8: Angular dependence of the particle fluxes f l20l ) for spin-0, spin-1/2 and spin-1 
fields for n = 1 and a* = 0.6. Here, 6 is the latitudinal angle, so that = 0, tt corresponds 
to the axis of rotation of the black hole and = 7r/2 is the equatorial plane. 

axis of the black hole, and cos^ = in the equatorial plane of the black hole. It can 
be seen that the angular distributions are quite different for the three cases. For the 
spin-0 case, we can again see the peaks in the particle flux. For small uj^ the angular 
distribution is uniform and spherically symmetric. However, for larger values of cj, the 
flux is concentrated in the equatorial plane of the black hole, as might be expected due 
to the rotation of the black hole. For the spin-1 case, the flux is again concentrated in 
the equatorial plane for large cj, but there is a more signiflcant concentration of emission 
along the axis of the black hole, for smaller values of uj. This difference in behaviour is due 
to the coupling of the rotation of the black hole with the electromagnetic fleld (which has 
non-zero spin). For the spin-1/2 case, the energy is again concentrated in the equatorial 
plane for large uj^ and there is a very slight spin-orbit coupling effect for small values of 

UJ. 
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Figure 9: Power fluxes f fT6l ) for spin-O, spin-1/2 and spin-l flelds as a function of urh^ for 
fixed n = 1 and varying a*. 




Figure 10: Power fluxes f fTGl ) for spin-0, spin-1/2 and spin-l flelds as a function of ujTh^ 
for flxed a* = 1 and varying n. 
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Next, we consider the power flux f fT6l ). and plot in flgures [9l and [TOl the power fluxes 
as functions of cjr/^, for flxed n and flxed a>^, respectively. Many of the features follow 
from those of the particle flux, as outlined above. In particular, emission at higher values 
of uj becomes more important as either n or a* increase. This is enhanced compared with 
the particle flux because the power in each mode involves multiplying the particle flux in 
that mode by uj. 

The angular dependence of the power flux f l2Tl ) is shown in flgure [H] for n = 1 and 
a* = 0.6. The spin-orbit coupling is particularly clear at small values of cj, and again the 




Figure 11: Angular dependence of the energy fluxes ( 12T1 ) for spin-0, spin- 1/2 and spin-1 
fields for n = 1 and = 0.6. Here, 9 is the latitudinal angle, so that ^ = 0, tt corresponds 
to the axis of rotation of the black hole and ^ = 7r/2 is the equatorial plane. 

emission is concentrated along the equatorial plane for large values of u. 

Finally, we consider the rate at which the black hole loses angular momentum. The 
fiux of angular momentum f fTTI l is shown in figures [l2] and [I3] for n = 1 and = 1 
respectively. The main features are as would be anticipated from our previous analysis of 
the particle and energy fiuxes. The angular dependence of the fiux of angular momentum 
has no experimental significance so we do not consider it here. 
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Figure 12: Angular momentum fluxes (fTTI) for spin-0, spin-1/2 and spin-1 as a function 
of UTh^ for flxed n = 1 and varying a*. 



spin spin 1/2 




Figure 13: Angular momentum fluxes (fTTI) for spin-0, spin-1/2 and spin-1 flelds as a 
function of ujTh^ for flxed a* = 1 and varying n. 
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5 Conclusions 



In this brief review we have considered the Hawking radiation of rotating brane black 
holes, as may be produced at the LHC. In references El ^ we have studied in detail 
the fluxes of particles, energy and angular momentum for particles of spin-0, spin-1/2 and 
spin-1, emitted on the brane. This is suflicient to compute the emission of all standard 
model particles on the brane. Our results will be of use for simulations of black hole 
events at the LHC. There are currently two such simulation codes, CHARYBDIS fL2\ [28] 
and CATFISH (see pO] for a recent comparison of these). 

To date both simulators only use detailed information about the 'Schwarzschild' 
phase of the black hole's evolution. It is hoped that our detailed results on the 'spin- 
down' phase will be incorporated into these simulators. This will enable the effects of the 
rotation of the black hole on LHC events to be accurately modelled. However, the analysis 
of the 'spin-down' phase is not yet complete: information about graviton emission (which 
will contribute 'missing energy' to LHC events [3l]) is required. We hope to return to 
this question in the near future. 
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